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Abstract 

c . 

It is believed that magnetic field lines are twisted and bend by turbulent motions in 
the Galaxy. Therefore, both Galactic synchrotron emission and thermal emission from 
dust reflects statistics of Galactic turbulence. Our simple model of Galactic turbulence, 
motivated by results of our simulations, predicts that Galactic disk and halo exhibit 
different angular power spectra. We show that observed angular spectra of synchrotron 
emission are compatible with our model. We also show that our model is compatible 
with the angular spectra of star-light polarization for the Galactic disk. Finally, we 
discuss how one can estimate polarized microwave emission from dust in the Galactic 
^ ' halo using star-light polarimetry. 
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Oh 1 Introduction 

6 

Exciting possibilities of measuring spectrum of polarized cosmic microwave background 
(CMB) fluctuations renewed the interest to the polarized fluctuations of Galactic origin. 
Synchrotron and dust are known to be the most important sources of polarized foreground 
radiation. Measurements of angular power spectra of such foregrounds are of great interest. 
So far, a large number of observations are available for angular power spectra of synchrotron 
emission and synchrotron polarization for Galactic disk and halo (see papers in de Oliveira- 
Costa & Tegmark 1999; see also references listed in §2 and papers in this volume). Those 
measurements revealed a range of power-laws, the origin of which has been addressed by 
Chepurnov (2002) and Cho & Lazarian (2002a; henceforth CL02) . The latter also addressed 
the origin of the observed angular spectrum of starlight polarization (Fosalba et al. 2002; 
henceforth FLPT). 

Interstellar medium is turbulent and Kolmogorov-type spectra were reported on the 
scales from several AU to several kpc (see Armstrong, Rickett, & Spangler 1995; Lazar- 
ian Sz Pogosyan 2000; Stanimirovic & Lazarian 2001; Lazarian & Esquivel 2003). Therefore 
it is natural to think of the turbulence as the origin of the fluctuations of the diffuse fore- 
ground radiation. Interstellar medium is magnetized with magnetic field making turbulence 
anisotropic. It may be argued that although the spectrum of MHD turbulence exhibits scale- 
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dependent anisotropy if studied in the system of reference denned by the local magnetic field 1 
(Goldreich & Sridhar 1995; Lithwick & Goldreich 2001; Cho & Lazarian 2002b,2003ab), in 
the observer's system of reference the spectrum will show only moderate scale-independent 
anisotropy. Thus from the observer's perspective Kolmogorov's description of interstellar 
turbulence statistics is acceptable in spite of the fact that magnetic field is dynamically im- 
portant and even dominant (see discussion in Lazarian & Pogosyan 2000; Cho, Lazarian, & 
Vishniac 2003). 

In this paper, we claim that MHD turbulence can explain observed angular spectra 
of Galactic synchrotron emission and starlight polarization. For this purpose, we use an 
analytical insight obtained in Lazarian (1992, 1995ab) and numerical results obtained in 
CL02.. We also discuss how we can estimate polarized microwave dust emission using star- 
light polarimetry. This problem is of great importance in view of recent interest to the 
foregrounds to the CMB polarization (see also review by Lazarian in this volume). 

2 Galactic synchrotron emission 
2.1 Summary of observations 

It is customary for CMB studies to expand the foreground intensity over spherical harmonics 
Y lm , 1(6, <p) = Y,i, m aimYim{0, <p) , and write the spectrum in terms of Ci = Em=-/ \ a im\ 2 / (2/ + 
1). The measurements indicate that angular power spectrum (C{) of the Galactic emission 
follows power law (C\ oc l~ a ) (see FLPT and references in §3). The multipole moment I is 
related to the angular scale 6 on the sky as I ~ ^jQradum ^ Qr ^ ^ 180° /6°). When the angular 
size of the observed sky patch (Ad x A9 in radian) is small, C\ is approximately the 'energy' 
spectrum of fluctuations (Bond & Efstathiou 1987; Hobson & Majueijo 1996; Seljak 1997). 
It is this power spectrum expressed in terms of wavenumber k (~ l Ad/ii) that is usually 
dealt with in studies of astrophysical turbulence (e.g. Stanimirovic & Lazarian 2001). 

Recent statistical studies of total synchrotron intensity include Haslam all-sky map at 
408 MHz (Haslam et al. 1982) that shows that the angular power index a (Ci oc l~ a ) is 
in the range between 2.5 and 3 (Tegmark & Efstathiou 1996; Bouchet, Gispert, & Puget 
1996). Parkes southern Galactic plane survey (Duncan et al. 1997) at 2.4 GHz suggests 
shallower slope: Giardino et al. (2002) obtained a ~ 2.05 after point source removal and 
Baccigalupi et al. (2001) obtained a ~ 0.8 to 2. On the other hand, Tucci et al. (2000) 
obtained a ~ 1.4 to 1.8 and Bruscoli et al. (2002) obtained a ~ 0.4 to 2.2 for the Galactic 
disk. Using Rhodes/HartRAO data at 2326 MHz (Jonas, Baart, & Nicolson 1998), Giardino 
et al. (2001a) obtained a ~ 2.43 for all-sky data and a ~ 2.92 for high Galactic latitude 
regions with \b\ > 20°. Giardino et al. (2001b) obtained a ~ 3.15 for high Galactic latitude 
regions with |6| > 20° from Reich & Reich (1986) survey at 1420 MHz. The rough tendency 
that follow from these data is that a which is close to 2 for the Galactic plane gets steeper 
(to ~ 3) for higher latitudes (Fig. |T^i). 

lr The anisotropy is present for AlfVenic and slow modes, while fast modes exhibit isotropy (see Cho & 
Lazarian 2002b). 
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Can we explain this tendency? In the next section, we construct a simple model that can 
explain the rough tendency of observed angular spectra. 



2.2 Basic idea of our simple model 

2.2.1 Ci for small angle limit 

In this section we show that, when the angle between the lines of sight is small (i.e. 9 < 
L/d max ), the angular spectrum Ci has the same slope as the 3-dimensional energy spectrum 
of turbulence. Here L is the typical size of the largest energy containing eddies, which is 
sometime called as outer scale of turbulence or energy injection scale, and d max is the distance 
to the farthest eddies (see Fig. Hb). 

To illustrate the problem consider the observations with lines of sight being parallel. The 
observed intensity is the intensity summed along the line of sight, r z . 

hD{r x ,r y ) = Jdr z i 3D (r) (1) 
= J dr z J dk x dk y dk z z 3D (k) e ikr . (2) 

Rearranging the order of summation and using / dr z e tkzTz = 8(k z ), we get 

h D (r x , r y ) = J dk x dk y i 3D (k x , k y , 0) e ik * r * +ik » r " , (3) 

which means Fourier transform of l2D{ r x, r y ) is i^D was mentioned earlier 

for small patches of sky C\ ~ \izoikx-, k y , 0)| 2 with I ~ k(n/A9) and k= (k 2 x + ky) 1/2 . 

The analysis of the geometry of crossing lines of sight is more involved, but for power-law 
statistics it follows from Lazarian & Shutenkov (1990) that if \izD\ 2 oc k~ m , then the 'energy' 
spectrum of l2D( r x, r y ) is also k~ m . Therefore, we have 

C,oc |i 3 D(^A,0)| 2 «^ m (4) 
in the small 9 limit. For Kolmogorov turbulence (I^dI 2 oc A;~ n / 3 ), we expect 

Q oc r 11/3 , if 9 < L/d max . (5) 

Note that / ~ n/9. 

2.2.2 Ci for large angle limit 

Following Lazarian & Shutenkov (1990), we can show that the angular correlation function 
for 9 > L/ d max is given by 



K{6) = J jdr 1 dr 2 }C(\v 1 -r 2 \), 



dr r/C(r) / dip ~ — - — , (6) 



sin 9 Jo Jo 9 
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Figure 1: (a) Comparison with observations. Widths of rectangles correspond to the Galactic 
latitude of the observed patches. Heights of rectangles denote either error bars or scatter of 
the slopes. ( Gia&8=Giardino et al. 2002; 5n/02=Bruscoli et al. 2002; etc.). (b) Upper plot: 
Small 9 limit (8 < L/d max ). The fluctuations along the entire length of the lines of sight 
are correlated. Lower plot: Large 6 limit [9 > L/d max ). Only points close to the observer 
are correlated, (c) Spectra for the disk (d max = 10 kpc; turbulence size L = 100 pc) and for 
the halo (d max = 1 kpc; turbulence size L = 100 pc). (d) Angular spectrum for the Galactic 
disk. We considered the geometry of the Galactic disk, (e) Angular spectra of degree of 
polarization. Fixed Distance means all stars are at the same distance of 10 kpc. Random 
Distance means random sampling of stars according to equation (JHJ). Zero points of C{6) and 
CVs are arbitrary, (f) Angular correlation function for polarized intensity at 1mm. From eq. 
(Unj). Figure (a)-(e) are from CL02. 

where we change variables from (r 1; r 2 ) to (r,ip), which is clear from Fig. lb. We accounted 
for the fact that the Jacobian of the transformation is rj sin#. We can understand 1/6 
behavior qualitatively from Fig. la. When the angle is large, the points along of the lines of 
sight near the observer are still correlated. These points extend from the observer over the 
distance oc 1/ sin (0/2). 

In the limit of 9 <C 1 we get the angular power spectrum C\ using Fourier transform: 

d ~ J J K(6)e- ll - e d8 x d8 y 
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I d9 6J (W)K(6) oc r\ (7) 



where 9 — (9 2 + 9 2 ) 1 / 2 , Jo is the Bessel function, and we use K{9) oc 9 1 . 

In summary, for Kolmogorov turbulence, we expect from equations (j5J) and (J7|) that 

f r 11 / 3 if / > I 

o <«{!- VKt (8) 

which means that the power index a of C\ is 2 — 1 < a < —11/3. The critical multipole 
moment l cr (~ 3d max /L) depends on the size of the large turbulent eddies and on the 
direction of the observation. If in the naive model we assume that turbulence is homogeneous 
along the lines of sight and that L ~ 100 pc, corresponding to a typical size of the supernova 
remnant, we get l cr ~ 300 for the Galactic disk with d max ~ 10 kpc. For the synchrotron 
halo, d 

max ~ 1 kpc (see Smoot 1999) and we get l cr ~ 30. 
2.3 Numerical results 

Can this model really explain observed behavior of synchrotron angular spectrum? To ad- 
dress the problem, we perform simple numerical calculations for the Galactic disk and halo. 
We obtain Ci using the relation 

K(cos9) = ^]T(2/ + l)QP,(cos#), (9) 

Q = ^J P l (cos9)K(cos9) d{cos9). (10) 

We use Gauss-Legendre quadrature integration (see Szapudi et al. 2001 for its application 
to CMB) to obtain C\. We numerically calculate the angular correlation function from 

K{cos9) = Jdn J dr 2 K{\r x -Y 2 \), (11) 

where |ri — r 2 | = r\ + r\ — 2r x r 2 cos# and assume that the spatial correlation function /C(r) 
follow Kolmogorov statistics: 

K(r)ocj . fr>L; (12) 

where /Co ~ L 2 ^ 3 is a constant. Fig. [I)p illustrates the agreement of our calculations with the 
theoretical expectations within the naive model of the disk and the halo from the previous 
section. 



2 Note that point sources would result in a ~ 0. 
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2.4 A more realistic model for Galactic disk 



The naive model discussed in the previous section predicts Z _1 spectrum for the Galactic disk 
for I < 300 (see Fig. Ht and the discussion below eq. (jHJ). However, as we discussed earlier, 
observations show angular spectra close to l~ 2 . In other words 2 < a < 3 which differs from 
naive expectations given by equation (JHJ). 

To make the spectrum closer to observations we need to consider more realistic models. 
First, synchrotron emission is stronger in spiral arms, and therefore we have more syn- 
chrotron emission coming from the nearby regions. Second, if synchrotron disk component is 
sufficiently thin, then lines of sight are not equivalent and effectively nearby disk component 
contributes more. Indeed, when we observe regions with low Galactic latitude, the effective 
line of sight varies with Galactic latitude. 

Suppose we observe a region with b = 10°. Then emission from d = 10 kpc is sub- 
stantially weaker than that from d = 100 pc, because the region with d — 10 kpc is 
d sin 10° ~ 1.7 kpc above the Galactic plane and, therefore, has weak emission. To in- 
corporate the effect of finite thickness (~ 100 pc) of the disk, we use the weighting function 
w(r) = 100pc/max(100pc, r sin 10°), which gives more weight to closer distances. The result- 
ing angular power spectrum (Fig. ^1) shows a slope similar to -2. 

For the halo, the simple model predicts that Ci oc Z™ 11 / 3 for I > 30, but observations 
provide somewhat less steep spectra. Is this discrepancy very significant? The spectrum of 
magnetic field is expected to be shallower than k~ n ^ 3 in the vicinity of the energy injection 
scale and at the vicinity of the magnetic equipartition scale. The observed spectrum also 
gets shallower if d max gets larger. For instance Beuermann, Kanbach, & Berkhuijsen (1985) 
reported the existence of thick radio halo that extends to more than ~ 3 kpc near the Sun. 
Finally, filamentary structures and point sources can make the spectrum shallower as well. 
Further research should establish the true reason for the discrepancy. 

3 MHD turbulence and Galactic Starlight Polarization 

Let us move to a different topic - starlight polarization. In this section, we will illustrate 
that MHD turbulence can explain the observed angular spectrum of starlight polarization. 

Polarized radiation from dust is an important component of Galactic foreground that 
strongly interferes with intended CMB polarization measurements (see Lazarian & Prunet 
2001). FLPT attempted to predict the spectrum of the polarized foreground from dust 
and obtained C\ ~ I for starlight polarization degree. They used polarization data from 
~5500 stars. The sample is strongly biased toward nearby stars (d < 2 kpc) within the 
Galactic plane. This spectrum is different from those discussed in the previous sections. To 
relate this spectrum to the underlying turbulence we should account for the following facts: 
a) the observations are done for the disk component of the Galaxy, b) the sampled stars are 
at different distances from the observer with most of the stars close-by. 

To deal with this problem we use numerical simulations again. We first generate a three 
(i.e. x,y, and z) components of magnetic field on a two-dimensional plane (4096 x 4096 grid 
points representing 10 kpc x 10 kpc), using the following Kolmogorov three-dimensional 
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spectrum: E^,£,[k) oc k~ u ^ 3 if k > ko, where ko ~ l/100pc. Our results show that the way 
how we continuously extend the spectrum for k < k does not change our results. 

To simulate the actual distribution of stars within the sample used in FLPT, we scatter 
our emission sources using the following probability distribution function: P(r) oc e ~ r / 1 - 5k P c _ 
The starlight polarization is due to the difference in absorption cross section of non-spherical 
grains aligned with their longer axes perpendicular to magnetic field (see a review by Lazarian 
2003). In numerical calculations we approximate the actual turbulent magnetic field by 
a superposition of the slabs with locally uniform magnetic field in each slap and assume 
that the difference in grain absorption parallel and perpendicular to magnetic field results 
in the 10% difference in the optical depths t\\ and r_i_ for a slab. We calculate evolution 
of Stokes parameters of the starlight within the slab and use the standard transforms of 
Stokes parameters from one slab to another (Dolginov, Gnedin, & Silantev 1996; see similar 
expressions in Martin 1972). 

We show the result in Fig.^. For comparison, we also calculate the degree of polarization 
assuming all stars are at the same distance of ~ 10 kpc. The result shows that, for a mixture 
of nearby and faraway stars, the slope steepens and gets very close to the observed one, i.e. 
-1.5. 



4 Estimation of Polarized Emission from Dust 

One of the possible ways to estimate the polarized radiation from dust at the microwave range 
is to measure star-light polarization and use the standard formulae relating polarization at 
different wavelength. The technique can be traced back to FLPT. In this section, we estimate 
polarized diffuse emission by dust in the high Galactic latitude halo. For our practical data 
handling, see Hildebrand et al. (1999). We use the data of optical depth (in fact, E(B — V)) 
and the degree of polarization by absorption P a bs,opticai for stars in the halo (provided by 
Terry Jones). 

When the optical depth is small, we have the following relation (see, for example, Hilde- 
brand et al. 2000): 

Pern, optical ~ P abs, optical / ^~ ■> (1*^) 

where P em ,opticai is the degree of polarization by emission and r is the optical depth (at optical 
wavelengths). We obtain polarization by emission at 1mm (P e m,m,m) using the relation 



Pem,mm Pem,optical 



Crn.ax. Cri 



n — n 

^max ^mm 



I 



n — n 

^max ^min 

~C _ r . 

^max ^mm 



(14) 

optical 



where C's are cross sections that depend on the geometrical shape (see, for example, the 
discussion in Hildebrand et al. 1999; see also Draine & Lee 1984) and dielectric function e 
(see Draine 1985) of grains. We assume grains are 2:1 oblate spheroids. 

We can obtain the emission at 1mm (I mm ) from that at 100 /im (hoo^m)'- 



'■ lOO^JTi 



(1mm / 100 fxm)~ p , (15) 
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where we use (3 ~ 1.7. For now, we assume that /loo^m ~ 0.85 MJy/sr (Boulanger & Perault 
1988; see also Draine & Lazarian 1998). From this, we obtain I mm « 0.017 MJy/sr. In the 
future, we can use an actual lOOfim map to get I mm . 

From P e m,mm and J mm , we can easily obtain the polarized intensity at 1mm (I po i, mm ): 

Ipol,mm Pem,mmlmm- (16) 

Fig. [Tf, obtained this way, shows I po i,mm- 

Last step will be the estimation of the angular spectrum C\. From CL02 (see also ^2.2j) . 
we obtain 

/A/- 1 if Z < 30 

\ 30 8 /MZ- u /3 if / > 30, U0 

where A is a constant that is determined by + l)Cj ~ 10~ 5 , where 10~ 5 is the rms 
fluctuation (see y-axis of Fig. Hf). Calculation yields A = 10~ 5 /30 2 . Therefore, 

r ~ J 10_7rl (MJy 2 /^ 2 ) if Z < 30 

lO- 4 /- 11 ^ (MJy 2 / S r 2 ) if Z > 30. [ ^ } 

When we prefer (fiK) to (MJy/sr), we can use the conversion factor c* ~ (l/x 2 )(10 4 fiK/ (MJy/sr)), 
where x ~ Z//56.8 GHz (Tegmark et al. 2000). For 1mm, x ~ 5 and c* ~ 400. Qc 2 gives 

f 10- 2 /- 1 (/i-ft' 2 ) if Z < 30 

\ 10/- n /3 (^K 2 ) if Z > 30. Uyj 

The particular normalization coefficients and the spectral index may vary from one region 
to another (see FLPT; CL02). 

We can summarize the procedure as follows: 



abs, optical * Pem,optical ^ Pem,mm 
1 100 fim * 1 mm 



d (at 1 mm) (20) 



5 Summary 

In this paper we have addressed the origin of spatial fluctuations of Galactic diffuse emission. 
We have shown that MHD turbulence with Kolmogorov spectrum can qualitatively explain 
the observed properties of synchrotron emission and starlight polarization. The variety of 
measured spatial spectra of synchrotron emission can be accounted for by the inhomogeneous 
distribution of emissivity along the line of sight arising from the structure of the Galactic disk 
and halo. Similarly, MHD turbulence plus inhomogeneous distribution of stars can explain 
the reported scaling of starlight polarization statistics. Although these complications do not 
allow to predict the exact slope of the measured fluctuations, our interpretation allows a 
valuable qualitative insight on what sort of change is to be expected. 

We have shown that one can estimate polarized diffuse micro-wave emission by dust 
using star-light. Our preliminary result shows that ST ~ [1(1 + l)Ci] l l 2 ~ 0(1) ixK for high 
Galactic latitude regions. 
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Evidently more systematic studies are required. Those studies will not only give insight 
into how to separate CMB from foregrounds, but also would provide valuable information 
on the structure of interstellar medium and the sources/energy injection scales of interstellar 
turbulence. (See also review by Lazarian, this volume.) 
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